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In this paper, Gaussianity of eigenmodes and non-Gaussianity in the Cosmic Microwave Back- 
ground (CMB) temperature fluctuations in two smallest compact hyperbolic (CH) models are in- 
vestigated. First, it is numerically found that the expansion coefficients of low-lying eigenmodes 
on the two CH manifolds behave as if they are Gaussian random numbers at almost all the places. 
Next, non-Gaussianity of the temperature fluctuations in the (l,m) space in these models is stud- 
ied. Assuming that the initial fluctuations are Gaussian, the real expansion coefficients bi m of the 
temperature fluctuations in the sky are found to be distinctively non-Gaussian. In particular, the 
cosmic variances are found to be much larger than that for Gaussian models. On the other hand, 
the anisotropic structure is vastly erased if one averages the fluctuations at a number of different 
observing points because of the Gaussian pseudo-randomness of the eigenmodes. Thus the domi- 
nant contribution to the two-point correlation functions comes from the isotropic terms described 
by the angular power spectra Ci. Finally, topological quantities: the total length and the genus of 
„Cj ■ isotemperature contours are investigated. The variances of total length and genus at high and low 

threshold levels are found to be considerably larger than that of Gaussian models while the means 
almost agree with them. 

■ PACS Numbers : 98.70.Vc, 98.80.Hw 
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In recent years, locally isortopic and homogeneous Friedmann-Robertson- Walker (FRW) models with non-trivial 
topology have attracted much attention. In the standard scenario, simply-connectivity of the spatial hypersurface is 
assumed for simplicity. However, the Einstein equations, being local equations, do not fix the global topology of the 
spacetime. In other words, a wide variety of topologically distinct spacetimes with the same local geometry described 
O ,■ by a local metric element remain unspecified (see [[y for review on the cosmological topology) . The determination of 
the global topology of the universe is one of the most important problem of the modern observational cosmology. 

For flat models without the cosmological constant, severest constraints have been obtained by using the COBE 
DMR data. The suppression of the fluctuations on scales beyond the topological identification scale L leads to the 
decrease of the angular power spectra Ci of the Cosmic Microwave Background (CMB) temperature fluctuations on 
large angular scales which puts a lower bound L>2400 h~ 1 Mpc (with h = Ho/100 kms _1 Mpc _1 ) for a compact flat 
. £h ] 3-torus model without the cosmological constant H[| . Similar constraints have been obtained for other compact flat 
^ i models {§. The maximum expected number of copies of the fundamental domain (cell) inside the last scattering 
surface is approximately 8 for the 3-torus model. 

In contrast, for low density models, the constraint could be considerably milder than the locally isotropic and 
homogeneous flat (Einstein-de-Sitter) models since a bulk of large-angle CMB fluctuations can be produced by the 
so-called (late) integrated Sachs- Wolfe effect (ISW) |||| which is the gravitational blueshift effect of the free streaming 
photons by the decay of the gravitational potential. As the gravitational potential decays in either A-dominant epoch 
or curvature dominant epoch, the free streaming photons with large wavelength (the light travel time across the 
wavelength is greater than or comparable to the decay time) that climbed a potential well at the last scattering 
experience blueshifts due to the contraction of the comoving space along the trajectories of the photons. Because the 
angular sizes of the fluctuations produced at late time are large, the suppression of the fluctuations on scale larger 
than the topological identification scale does not lead to a significant suppression of the large-angle power if the ISW 
effect is dominant. Recent works J?]-^0| have shown that the large-angle power (2 < Z < 20) are completely consistent 
with the COBE DMR data for compact hyperbolic (CH) models which include a small CH orbifold, the Weeks and 
the Thurston manifolds with volume 0.72, 0.94 and 0.98 in unit of the cube of the curvature radius, respectively. 
Note that the Weeks manifold is the smallest and the Thurston manifolds is the second smallest in the known CH 
manifolds. For instance, the number of copies of the fundamental domain inside the last scattering surface at present 
is approximately 190 for a Weeks model with f2o = 0.3. 

If the space is negatively curved, for a fixed number of the copies of the fundamental domain inside the present 
horizon, the large-angle fluctuations can be produced much effectively. In negatively curved spaces (hyperbolic spaces), 
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trajectories of photons subtend a much smaller angle in the sky for a given scale. In other words, for a given angle 
of a pair of two photon trajectories, the physical distance of the trajectories is much greater than that in a flat 
space. Therefore, even if there is a number of copies of the fundamental domain which intersect the last scattering 
surface, the number of copies which intersect the wave front (a sphere with z = const.) of the free streaming photons 
is exponentially decreased at late time when the large- angle fluctuations are produced due to the ISW effect. 

However, one may not be satisfied with the constraints using only the angular power spectrum C; since it contains 
only isotropic information of the ensemble averaged temperature fluctuations JH| . If they have anisotropic structures, 
non-Gaussian signatures must be revealed. In fact, the global isotropy of the locally isotropic and homogeneous 
FRW models is generally broken. For instance, a flat 3-torus obtained by identifying the opposite faces of a cube is 
obviously anisotropic at any points. Thus the temperature fluctuations averaged over the initial conditions in these 
multiply-connected FRW models are no longer SO(3) invariant at a certain point. The temperature fluctuations on 
the sky are written in terms of (real) spherical harmonics Q; m (n) as 

AT - 

-jT( n )=Zj b l™Qhn(n). (I) 
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If the distribution functions of the real expansion coefficients &; m are SO(3) invariant, the temperature fluctuations 
must be Gaussian provided that &z m 's are independent random numbers [ ft2[ . Therefore, the temperature fluctuations 
at a certain point in the multiply-connected FRW models are not Gaussian if 6; m 's are independent. 

For the simplest flat 3-torus models (without rotations in the identification maps) which are globally homogeneous, 
it is sufficient to choose one observing point and estimate how the power is distributed among the m's for a given 
angular scale I in order to see the effect of the global anisotropy. However, in general, one must consider an ensemble 
of fluctuations at different observing points because of the spatial (global) inhomogeneity. Previous analyses have not 
fully investigated the dependence of the temperature fluctuations on choice of the observing points. 

Lack of analytical results on the eigenmodes makes it difficult to investigate the nature of the temperature fluctua- 
tions in CH models. However, we may expect a high degree of complexity in the eigenmodes since the corresponding 
classical systems (geodesic flows) are strongly chaotic. In fact, it has been numerically found that the expansion 
coefficients of the low-lying eigenmodes on the Thurston manifold at the point where the injectivity radius is maximal 
are Gaussian pseudo-random numbers ]T^ | which supports the previous analysis of the excited states (higher modes) 
of a two-dimensional asymmetrical CH model [Q. We have put a prefix "pseudo" since the eigenmodes are actually 
constrained by the periodic boundary conditions. These results imply that the statistical properties of the eigenmodes 
on CH spaces (orbifolds and manifolds) can be described by random-matrix theory (RMT) |lj,[f(| . An investigation of 
the dependence of the property on the observing points is also important since CH spaces have symmetries (isometric 
groups) which may veil the random feature of the eigenmodes. In this paper, a detailed analysis on the statistical 
property of low-lying eigenmodes on the Weeks and the Thurston manifolds is conducted. 

Assuming that the eigenmodes are Gaussian, one can expect that the anisotropic structure in the (l,m) space is 
vastly erased when one averages the fluctuations over the space. This seems to be a paradox since the CH spaces 
are actually globally anisotropic. However, one should consider a spatial average of fluctuations with different initial 
conditions if one believes the Copernican principle that we are not in the center of the universe. Even if the space is 
anisotropic at a certain point, the averaged fluctuations may look isotropic by considering an ensemble of fluctuations 
at all the possible observing points. Note that the eigenmodes on CH spaces have no particular directions if they are 
Gaussian. 

If the initial fluctuations are constant for each eigenmode, as we shall see, the Gaussian randomness of the tem- 
perature fluctuations can be solely attributed to the Gaussian pseudo-randomness of the eigenmodes. In this case, 
the Gaussian randomness of the temperature fluctuations has its origin in the geometrical property of the space ( Ge- 
ometric Gaussianity). Choosing an observing point is equivalent to fixing a certain initial condition. However, it is 
much natural to assume that the initial fluctuations are also random Gaussian as the standard inflationary scenarios 
predict. Then the temperature fluctuations may not obey the Gaussian statistics because they are written in terms of 
products of two different independent Gaussian numbers rather than sums while they remain almost spatially isotropic 
if averaged over the space. 

In this paper, Gaussianity of eigenmodes and non-Gaussianity in the CMB for two smallest CH models (the Weeks 
and the Thurston models) are investigated. In Sec. II, numerical results on Gaussianity of eigenmodes are shown and 
we discuss to what extent the results are generic. In Sec. Ill, we study the non-Gaussian behavior of the temper- 
ature fluctuations in the (l,m) space. In Sec. IV, topological quantities (total length and genus) of isotemperature 
contours are numerically simulated for studying the non-Gaussian behavior in the real space. Finally, we summarize 
our conclusions in Sec. V. 
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II. GEOMETRIC GAUSSIANITY 



In locally isotropic and homogeneous FRW background spaces, each type (scalar, vector and tensor) of first-order 
perturbations can be decomposed into a decoupled set of equations. In order to solve the decomposed linearly 
perturbed Einstein equations, it is useful to expand the perturbations in terms of cigcnmodes of the Laplacian which 
satisfies the Helmholtz equation with certain boundary conditions, 

(V 2 + k 2 )u k (x) =0, (2) 

since each eigenmode evolves independently in the linear approximation. Then one can easily see that the time evo- 
lution of the perturbations in the multiply-connected locally isotropic and homogeneous FRW spaces coincide with 
that in the FRW spaces while the global structure of the background space is described solely by these eigenmodes. 

Unfortunately, no analytical expressions of eigenmodes on CH spaces have been known. Nevertheless, the corre- 
spondence between classical and quantum mechanics may provide us a clue for understanding the generic property 
of the eigenmodes. If one recognizes the Laplacian as the Hamiltonian in a quantum system, each eigenmode can be 
interpreted as a wavefunction in a stationary state. Because classical dynamical systems (^geodesic flows) on CH 
spaces are strongly chaotic (or more precisely they are K-systems with ergodicity, mixing and Bernoulli properties 
), one can expect a high degree of complexity for each eigenstate. Imprint of the chaos in the classical systems may 
be hidden in the quantum counterparts. In fact, in many cases, the short-range correlations observed in the eigen- 
values (energy states) have been found to be consistent with the universal prediction of RMT for three universality 
classes:the Gaussian orthogonal ensemble(GOE), the Gaussian unitary ensemble(GUE) and the Gaussian symplectic 
ensemble (GSE) fj|||J|]. In our case the statistical properties are described by GOE (which consist of real symmetric 
NxN matrices H which obey the Gaussian distribution cx cxp (— TriJ 2 / (4a 2 )) (where a is a constant) as the systems 
possess a time-reversal symmetry. RMT also predicts that the squared expansion coefficients of an eigenstate with re- 
spect to a generic basis are distributed as Gaussian random numbers []l8| . Unfortunately, no analytic forms of generic 
bases(=eigenmodes) are known for CH spaces which seems to be an intractable problem. However, if the eigenmodes 
are continued onto the universal covering space by the periodic boundary conditions, they can be written in terms of 
a "generic" basis on the universal covering space (=3-hyperboloid H 3 ). In pseudospherical coordinates (R,x,d,(j)), 
the eigenmodes are written in terms of complex expansion coefficients ^ v i m and eigenmodes on the universal covering 
space, 

U v = ^ £rslm Xvl (x)Y~l, n (0,4>), (3) 
Im 



where v = \Jk 2 — 1, X v i and Yi m denote the radial eigenfunction and (complex) spherical harmonic on the pseudo- 
sphere with radius R, respectively. Then the real expansion coefficients a„/ m are given by 

a^oo = — Im^oo), a»io — V^^^zo), 
a vlrn = v / 2Re(^ I/ j m ), m > 0, 

a v im = -V2Im(£„ ; _ m ), m < 0, (4) 

where 



c v i 



(l + Re(F(^,0))' 



v(i + + n-vi) 

F ^ l) = v(vi) r(i-* + ir (5) 



In this paper, the low-lying eigenmodes (k < 13) on the Weeks and Thurston manifolds are numerically computed by 
the direct boundary element method. The identification matrices of the Dirichlet domains are obtained by a computer 
program "SnapPea" by J. R. Weeks |19[j . The computed eigenvalues are well consistent with that in the previous 
literature [ p^| , p0[ . The estimated errors in k are within 0.01. However, the last digits in k may be incorrect. a„/ m 's 
can be promptly obtained after the normalization and orthogonalization of these eigenmodes. The orthogonalization 
is achieved at the level of 10 -3 to 10 -4 (for the inner product of the normalized eigenmodes) which implies that each 
eigenmode is computed with relatively high accuracy. In Fig.l and Fig. 2, one can see a high degree of complexity in 
the lowest eigenmodes on the Poincare ball which is isometric to the universal covering space H 3 whose coordinates 
are given by 
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Weeks manifold 




FIG. 1. The lowest eigenmode k = 5.268 on the Weeks manifold continued onto the Poincare ball and the boundaries of the 
copied Dirichlet domains (solid curves) plotted on a slice z~0. 
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FIG. 3. Plots of a„ !m 's which are ordered as 1(1 + 1) + m + 1, < I < 20 for eigenmodes fc = 5.268(left) and k = 12.789(right) 
on the Weeks manifold at a point which is randomly chosen. 



x = R tanh ~ sin 9 cos 0, y — R tanh ^ sin 9 sin 0, z = i?tanh ^ cos#. (6) 

Replacing tanh ^ by tanh% for each coordinate, one obtains the Klein (projective) coordinates. In the Poincare 
coordinates, angles of geodesies coincide with that of Euclidean ones. In the Klein coordinates, all geodesies are 
straight lines while angles does not coincide with that of Euclidean ones. In what follows R is normalized to 1 without 
loss of generality. 

In Fig. 3, one can see that the distribution of a y ; m 's which are ordered as 1(1 + 1) + m + 1 are qualitatively random. 
In order to estimate the randomness quantitatively, we consider a cumulative distribution of 

, \o>vlm Q'v\ / n \ 
Ovlm — n { I ) 

where a v is the mean of a y ; m 's and a 2 is the variance. If a y ; m 's are Gaussian then b v i m 's obey a \ 2 distribution 
P(x) — (l/2) 1 l 2 T(l/2)x~ 1 / 2 e'~ x / 2 with 1 degree of freedom. To test the goodness of fit between the the theoretical 
cumulative distribution I(x) and the empirical cumulative distribution function In(x), we use the Kolmogorov- 
Smirnov statistic Dn which is the least upper bound of all pointwise differences \In(x) — I(x)\ plfl, 



D N = aup\I N (x)-I(x)\. (8) 



In(x) is defined as 




x < yi, 

I N (x) = { j/N, y 3 <x< y j+1 , j = 1, 2, . . . , N-l, (9) 

Vn < x, 

where yx < < ■ ■ ■ < yN are the computed values of a random sample which consists of N elements. For random 
variables Dn for any z > 0, it can be shown that the probability of <d is given by |^2| 

lim P(D N <d = zN~ 1/2 ) = L(z), (10) 

N— »oo 

where 

oo 

i(z) = l-2^(-ir 1 e- 2 ^ 2 . (11) 

3=1 

From the observed maximum difference Djy = d, we obtain the significance level ao = I — P which is equal to the 
probability of > d. If au is found to be large enough, the hypothesis In(x)^I(x) is not verified. The significance 
levels ajv for 0< I < 20 for eigenmodes k < 13 on the Thurston manifold are shown in table 1. The agreement with 
the RMT prediction is fairly good for most of eigenmodes which is consistent with the previous computation in |l3| . 
However, for five degenerated modes, the non-Gaussian signatures are prominent (in Jl3[ , two modes in (k< 10) have 
been missed). Where does this non-Gaussianity come from? 

First of all, we must pay attention to the fact that the expansion coefficients a v \ m depend on the observing point. 
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TABLE I. Eigenvalues k and the corresponding significance levels od for the test of the hypothesis In{x) 7^ I(x) for the 
Thurston manifold. The injectivity radius is maximal at the base point. 



In mathematical literature the point is called the base point. For a given base point, it is possible to construct a 
particular class of fundamental domain called the Dirichlet (fundamental) domain which is a convex polyhedron. A 
Dirichlet domain fl(x) centered at a base point x is defined as 

Sl(x) = p| H(g, x) , H(g, x) = {z\d(z, x) < d(g(z), x)}, (12) 
9 

where g is an element of a Klcinian group T(a discrete isometry group of PSL(2, C)) and d(z, x) is the proper distance 
between z and x. 

The shape of the Dirichlet domain depends on the base point but the volume is invariant. Although the base point 
can be chosen arbitrarily, it is a standard to choose a point Q where the injective radius [] is locally maximal. More 
intuitively, Q is a center where one can put a largest connected ball on the manifold. If one chooses other point as 
the base point, the nearest copy of the base point can be much nearer. The reason to choose Q as a base point is that 
one can expect the corresponding Dirichlet domain to have many symmetries at Q |p3f . 

As shown in Fig. 4, the Dirichlet domain at Q has a Z2 symmetry (invariant by 7r-rotation) if all the congruent 
faces are identified. Generally, congruent faces are distinguished but it is found that these five modes have exactly 
the same values of eigenmodes on these congruent faces. Then one can no longer consider a„/ m 's as "independent" 
random numbers. Choosing the invariant axis by the 7r-rotation as the z-axis, a v imS are zero for odd m's which leads 
to the observed non-Gaussian behavior. It should be noted that the observed Z2 symmetry is not the subgroup of 
the isometry group (or symmetry group in mathematical literature) D2 (dihedral group with order 2) of the Thurston 
manifold since the congruent faces must be actually distinguished in the manifold^] 

Thus the observed non-Gaussianity is caused by a particular choice of the base point. However, in general, the 
chance that we actually observe any symmetries (elements of the isometry group of the manifold or the finite sheeted 



1 The injective radius of a point p is equal to half the length of the shortest periodic geodesic at p. 

The observed Z2 symmetry is considered to be a "hidden symmetry" which is a symmetry of the finite sheeted cover of 
the manifold (which tessellates the manifold as well as the universal covering space). For instance, the Dirichlet domain of 
the Thurston manifold can be tessellated by four pieces with three neighboring kite-like quadrilateral faces and one equilateral 
triangle on the boundary and seven faces which contain the center as a vertex. By identifying the four pieces (by a tetrahedral 
symmetry), one obtains an orbifold which has a Z2 symmetry. 
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FIG. 4. A Dirichlet domain of the Thurston manifold in the Klein coordinates viewed from opposite directions at Q where 
the injectivity radius is locally maximal. The Dirichlet domain has a Z2 symmetry(invariant by 7r-rotation)at Q. 

cover of the manifold) is expected to be very low. Because a fixed point by an element of the isometric group is 
either a part of 1-dimensional line (for instance, an axis of a rotation) or an isolated point (for instance, a center of 
an antipodal map). 

In order to confirm that the chance is actually low, the KS statistics old of a„; m 's are computed at 300 base points 
which are randomly chosen. As shown in table 2, the averaged significance levels < old > are remarkably consistent 
with the Gaussian prediction, la of old are found to be 0.26 to 0.30. 

Next, we apply the run test for testing the randomness of a„; m 's where each set of a„/ m 's are ordered as l(l+l)+m+l 
(see Suppose that we have n observations of the random variable U which falls above the median and n 

observations of the random variable L which falls below the median. The combination of those variables into 2n 
observations placed in ascending order of magnitude yields 



Each underlined group which consists of successive values of U or L is called run. The total number of run is called 
the run number. The run test is useful because the run number always obeys the Gaussian statistics in the limit 
n — * oo regardless of the type of the distribution function of the random variables. As shown in table 3, averaged 
significance levels <a r > are very high (la is 0.25 to 0.31). Thus each set of a„; m 's ordered as 1(1 + 1) + m + 1 can be 
interpreted as a set of Gaussian pseudo-random numbers except for limited choices of the base point where one can 
observe symmetries of eigenmodes. 

Up to now, we have considered I and m as the index numbers of a u i m at a fixed base point. However, for a fixed 
(l,m), the statistical property of a set of a„/ m 's at a number of different base points is also important since the 
temperature fluctuations must be averaged all over the places for spatially inhomogeneous models. From Fig. 5, one 
can see the behavior of m-averaged significance levels 



which are calculated based on 300 realizations of the base points. It should be noted that each a v i m at a particular 
base point is now considered to be "one realization" whereas a choice of I and m is considered to be "one realization" 
in the previous analysis (table 1). The agreement with the RMT prediction is considerably good for components l> 1. 
For components 1 = 1, the disagreement occurs for only several modes. However, the non-Gaussian behavior is distinct 
in 1 = components. What is the reason of the non-Gaussian behavior for I = 0? Let us estimate the values of the 
expansion coefficients for 1 = 0. In general, the complex expansion coefficients ^ v i m can be written as, 
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(13) 




(14) 



For 1 = 0, the equation becomes 




(15) 
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TABLE II. Eigenvalues k and corresponding averaged significance levels < an > based on 300 realizations of the base points 
for the test of the hypothesis In{x) / I(x) for the Weeks and the Thurston manifolds. 
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9,881b 


0.51 


12.789 


0.48 


9.904 


0.52 


12.795 


0.50 


10.335a 


0.54 






9.984 


0.49 


12.806 


0.51 


10.335b 


0.51 






10.358 


0.53 


12.897a 


0.57 


10.452a 


0.53 






10.686a 


0.51 


12.897b 


0.55 



TABLE III. Eigenvalues k and corresponding averaged significance levels < a r > for the test of the hypothesis that the a„i m 's 
are not random numbers for the Weeks and Thurston manifolds. a r 's at 300 points which are randomly chosen are used for 
the computation. 
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FIG. 5. Plots of m-averaged significance levels qd(i' 1 I) based on 300 realizations for the Weeks and the Thurston manifolds 
(0<i<20 and k < 13). n denotes the index number which corresponds to an eigenmode Uk where the number of eigenmodes 
less than k is equal to n (k(n = 1) is the lowest non-zero eigenvalue). The accompanying palettes show the correspondence 
between the level of the grey and the value. 



Taking the limit xo ~ * 0, one obtains, 

2itu v (0)i 

&/00 = ■ (16) 

Thus a^oo can be written in terms of the value of the eigenmode at the base point. As shown in Fig.l, the lowest 
eigenmodes have only one "wave" on scale of the topological identification scale L (which will be defined later on) 
inside a single Dirichlet domain which implies that the random behavior within the domain may be not present. 
Therefore, for low-lying eigenmodes, one would generally expect non-Gaussianity in a set of a„oo 's. However, for 
high- lying eigenmodes, this may not be the case since these modes have a number of "waves" on scale of L and they 
may change their values locally in a almost random fashion. 

The above argument cannot be applicable to a v i m 's for I ^ where X„i approaches zero in the limit Xo —> while 
the integral term 

/ 'u„(xo,M)ijm(M)dn (17) 

also goes to zero because of the symmetric property of the spherical harmonics. Therefore o„j m 's cannot be written in 
terms of the local value of the eigenmode for I ^ 0. For these modes, it is better to consider the opposite limit xo —* 00 • 
It is numerically found that the sphere with very large radius xo intersects each copy of the Dirichlet domain almost 
randomly (the pulled back surface into a single Dirichlet domain chaotically fills up the domain) . Then the values of 
the eigenmodes on the sphere with very large radius vary in an almost random fashion. For large xoj we bave 

A^(xo)°ce- 2 *°+^<)\ (18) 

where <f>(y, I) describes the phase factor. Therefore, the order of the integrand in Eq. ( |l4| ) is approximately e~ 2x ° 
since Eq. (jlj) does not depend on the choice of xo- As the spherical harmonics do not have correlation with the 
eigenmode u v {xoi 0, </>), the integrand varies almost randomly for different choices of (l,m) or base points. Thus we 
conjecture that Gaussianity of a^m's have their origins in the chaotic property of the sphere with large radius in CH 
spaces. The property may be related to the classical chaos in geodesic fiowsR . 



3 If one considers a great circle on a sphere with large radius, the length of the circle is very long except for rare cases in 
which the circle "comes back" before it wraps around in the universal covering space. Because the long geodesies in CH spaces 
chaotically (with no particular direction and position) wrap through the manifold, it is natural to assume that the great circles 
also have this chaotic property. 
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So far we have seen the Gaussian pseudo-randomness of the a„j m 's. Let us now consider the statistical properties 
of the expansion coefficients. As the eigenmodes have oscillatory features, it is natural to expect that the averages are 
equal to zero. In fact, the averages of < a„; m > 's over < I < 20 and —l<m<—l and 300 realizations of base points for 
each i/-mode are numerically found to be 0.006 ± 0.04 - 0.02(1<t) for the Weeks manifold, and 0.003 ± 0.04 - 0.02(1ct) 
for the Thurston manifold. Let us next consider the ^-dependence (fc-dependence) of the variances Var{a v i rn ). In 
order to crudely estimate the ^-dependence, we need the angular size 59 of the characteristic length of the eigenmode 
u v at xo U| 

S6 2 a 16^ Vol(M) 

k 2 (smh(2(xo + r ave )) - sinh(2(x o ^ave ))-4r 

where Vol(M) denotes the volume of a manifold M and r ave is the averaged radius of the Dirichlet domain. There is 
an arbitrariness in the definition of r ave . Here we define r ave as the radius of a sphere with volume equivalent to the 
volume of the manifold, 

Vol(M) = 7r(sinh(2r al , e ) - 2r ave ), (20) 

which does not depend on the choice of the base point. The topological identification length L is defined as L — 2r ave . 
For the Weeks and the Thurston manifold, L — 1.19 and L — 1.20 respectively. From Eq. (|l9|), for large xot one can 
approximate u v (x ) ~ u' u {x'o) by choosing an appropriate radius x'o which satisfies v~ 2 exp(— 2x ) = v'~ 2 exp(— 2x' )- 
Averaging Eq. (n4|) over Vs and m's or the base points, for large Xoi one obtains, 

(i^^i 2 )~ exp i~f;| (i^i 2 ), (21) 

exp(-2x^) 

which gives (|^ w i m | 2 ) ~ v~ 2 . Thus the variance of a^imS is proportional to v~ 2 . The numerical results for the two 
CH manifolds shown in Fig. 6 clearly support the v~ 2 dependence of the variance. 

As we have seen, the property of eigenmodes on general CH manifolds is summarized in the following conjecture: 



Conjecture: Except for the base points which are too close to any fixed points by symmetries, for a fixed v, a set 
of the expansion coefficients a v i m over (l,m) 's can be considered as Gaussian pseudo-random numbers. For a fixed 
(vim) (I > 0), the expansion coefficients at different base points that are randomly chosen can also be considered as 
Gaussian pseudo-random numbers. In either case, the variance is proportional to v~ 2 and the average is zero. 



III. NON-GAUSSIANITY IN OBSERVABLE ANGULAR POWER SPECTRA 



As mentioned in the last section, perturbations in CH models are written in terms of linear combinations of 
eigenmodes and the time evolution of the perturbations. Because the time evolution of the perturbations coincides 
with that in open models, once the expansion coefficients ^ v i m (or a„i m ) are given, the evolution of perturbations in 
CH models can be readily obtained. 
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If one assumes that the perturbation is a adiabatic scalar type without anisotropic pressure, and the subhorizon 
effects such as acoustic oscillations of the temperature and the velocity of the bulk fluid, and the effect of the 
radiation contribution at high z are negligible, the time evolution of the growing mode of the Newtonian curvature $ 
is analytically given as (see e.g. p4| , ^5f ) 

, . 5(sinh 2 77 — 3t? sinh n + 4 cosh 77 — 4) . 

^ = ~ (C0sh,! 1)3 ^ > ^ 

where rj denotes the conformal time. In terms of <&, the temperature fluctuation in the sky are written as 

AT(n) ^ 

— ^ — = 2 , a '™^m(n) 

Im 

= ^^(0)^^(770)^(11), (23) 



vim 



where 



1 f Vo d<f> 

F„j(Tto)= -~$(<n*)X vl (r)Q-i 1 *)-2 dri—X u i(r) -v)- (24) 

Here $j,(0) is the initial value of the curvature perturbation and 7/» and 770 are the conformal time of the last scattering 
and the present conformal time, respectively. The angular power spectrum C\ is defined as 

1 

(2i + i)c,= <l a <™! 2 > 

tn=—l 

(|^™| 2 )|^(%)| 2 , (25) 



E 



v(v 2 + l)Vol(M) 



where V<s>{v) is the initial power spectrum. It should be noted that the above formula converges to that of open 
models in the short- wavelength limit (summation to integration) provided that < |^; m | 2 > is proportional to v~ 2 . 
The reason is as follows: Let us denote the number of eigenmodes with eigenvalues equal to or less than v by N(v). 
In the short-wavelength limit v» 1 one can use Weyl's asymptotic formula which leads to 



dN{v) Vol(M) 



dv 2tt 2 



v\ (26) 



Thus the v 2 dependence in Eq.(26) is exactly cancelled out by the v~ 2 dependence of eigenmodes. In what follows 
we assume the extended Harrison-Zcl'dovich spectrum, i.e. V^{v) = Const. (in the flat limit, it converges to the scale 
invariant Harrison-Zel'dovich spectrum) as the initial power spectrum. 

In estimating the temperature correlations, the non-diagonal terms (I/f or m!) may not be negligible if the 
background spatial hypersurface is not isotropic, in other words, the angular power spectrum C; may not be sufficient 
in describing the temperature correlations since Ci provides us with only an isotropic information of statistics of the 
correlations. However, this is not the case for CH models to which the conjecture proposed in Sec. II is applicable. 
Based on the Copernican principle, it is not likely that we are at the center of any symmetries. Therefore, in order 
to statistically estimate the temperature correlations in the globally inhomogeneous background space, one has to 
consider an ensemble of fluctuations with different initial conditions at different places (or base points) with different 
orientations. Almost all the anisotropic information is lost in the spatial averaging process since the eigenmodes are 
Gaussian. 

As shown in Fig. 7, for 300 realizations of observing points(left), the averaged absolute values of the off-diagonal 
elements in unit of diagonal elements are very small (~ 0.016) whereas their contributions seem to be not negligible 
(~ 0.25) at one particular observing point(right) where one can observe a symmetry of the Dirichlet domain. Thus 
the statistical property of the temperature correlation can be estimated by using C;'s provided that the eigenmodes 
are Gaussian which validates the previous analyses using Cz's for constraining the CH models @ 0. The spatial 
averaging process^ must be taken into account since there is no reason to believe that we are in the center of any 



4 In general, one should include an averaging process over different choices of orientation of coordinates as well as an averaging 
process over different choices of the observing point. Nevertheless, the Gaussian conjecture in Sec. II implies that the eigenmodes 
on CH spaces are "SO(3) invariant" if averaged all over the space. Therefore, omission of the averaging procedure for 
different orientations of coordinates make no difference. 
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FIG. 7. Contributions of non-diagonal terms in the temperature correlations in unit of diagonal terms which are defined 
as fjl^i = j < ai m a*ir m i > [ /y< |a; m | 2 >< |a;'m'| 2 > f° r the Thurston model with f2o=0.3. The four-dimensional space 
(I, m, I', m') is represented in the two-dimensional space as (n, n') = (1(1 + 1) + m + 1, Z'(Z' + 1) + m' + 1) for 2 < Z < 10, — I < m < I 
and 2 < Z' < 10, — Z' < m' < l ' . ffl^i's are represented by the level of grey shown in the accompanying palettes. The left figure 
represents fjl^i 's averaged over 300 realizations of the base points with infinite number of initial conditions for the Newtonian 
curvature. The right figure represents ' s a t a base point where the injective radius is maximal with infinite number of initial 
conditions. The computation is based on 36 eigenmodes(Zc < 13) that are numerically obtained by using the direct boundary 
element method. The averaged values of the non-diagonal fl^/'s (Z^Z' or m^m')a,re 0.016(left) and 0.25(right). 



symmetries. 

If the initial conditions satisfy (<lv(0))~ 2 oc v(v 2 + 1) that corresponds to the extended Harrison-ZePdovich spectrum 
, then Eq.(|2^) tells us that the temperature fluctuation is Gaussian since it is equal to a sum of Gaussian (pseudo- 
random numbers at almost all the observing points. In this case, the Gaussian randomness of the temperature 
fluctuations in CH models can be solely attributed to the geometrical property of the space (geometric Gaussianity) 
which may be related to the deterministic chaos of the corresponding classical system. In other words, the Gaussian 
randomness can be explained in terms of the classical physical quantities without considering the initial quantum 
fluctuations provided that the above conditions are initially (deterministically) satisfied. 

However, it is much natural to assume that 4> l/ (0)'s are also random Gaussian as in the inflationary scenarios in 
which Gaussianity (on large scales) of the temperature fluctuations has its origin in Gaussianity of the initial quantum 
fluctuations because the angular powers are generally similar to the extended Harisson-Zel'dovich spectrum. Then the 
statistical properties of the temperature fluctuations are determined by the sum of the products of the two independent 
Gaussian random numbers (the initial fluctuations and the expansion coefficients of the eigenmodes). 

Let us calculate the distribution function F(Z, az) of a product of two independent random numbers X and Y that 
obey the Gaussian (normal) distributions N(X;0,ax) and N(Y;0,ay), respectively. 



N(X- l x,a) = 



1 



-c 



-(X-^) 2 /2a 2 



2vrcr 



(27) 



Then F(Z — XY, az) is readily given by 



dY 



F(Z, a z )=2 N(Z/Y, 0, a x )N(Y, 0, a Y ) 

Jo Y 

_^_*,(J£L 

7T(jx <ry ^crxcy 



(28) 



where Kq(z) is the modified Bessel function. The average of Z is zero and the standard deviation satisfies az = axay- 
As is well known, Kq(z) is the Green function of the diffusion equation with sources distributed along an infinite line. 
Although Kq(z) is diverged at z — its integration over (—00,00) is convergent. From the asymptotic expansion of 
the modified Bessel function 



*o(*W S e- 



V 



l 2 -3 2 



l!8z 2\(Szf 



l 2 ■ 3 2 ■ 5 2 

3!(8z) 3 



z » 1, 



(29) 
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FIG. 8. On the left, the distribution function F(Z, 1) for a product of two random Gaussian numbers is plotted in solid 
curves. On the right, the distribution function G(Z, 1) (la — 1) of a sum of two random variables that obey F(Z, 1/V2). The 
dashed curves represent the Gaussian distribution N(Z; 0, 1). 



one obtains in the lowest order approximation, 

F(Z,,)~^i=e-l z ^ Z»l. (30) 

Thus F (Z, a) is slowly decreased than the Gaussian distribution function with the same variance in the large limit . 
One can see the two non-Gaussian features in Fig.8(left):the divergence at Z — > and the slow convergence to zero 
at Z — > oo. The slow convergence is an important feature, as we shall see, in distinguishing the non-Gaussian models 
with the Gaussian ones. In the modest region 0.4 < \Z\ < 2.4, F(Z,l) is much less than N(Z, Q, 1). Generally, 
the temperature fluctuation is written as a sum of the random variables Zi which obeys the distribution function 
F(Zi,<JZi) for a fixed set of cosmological parameters. For large-angle fluctuations, only the eigenmodes with large 
wavelength (=2-7r/fc)can contribute to the sum. Due to the finiteness of the space, the number of eigenmodes which 
dominantly contribute to the sum is finite. Therefore, the fluctuations are distinctively non-Gaussian. For small-angle 
fluctuations, the number of eigenmodes that contribute to the sum becomes so large that the distribution function 
converges to the Gaussian distribution as the central limit theorem implies. One can see from Fig. 8 (right) that the 
distribution function G(W, 1) of W — Zl + Z2 where both Zl and Z2 obey F(Z, \/2) is much similar to the Gaussian 
distribution N(Z, 0,1) than F(W, 1) in the modest region. 

Now let us see the non-Gaussian features of the observable angular power spectrum C; assuming that the initial 
fluctuations are Gaussian. First of all, we define a statistic x 2 = (21 + l)Ci/Ci where 

i 

(2f + l)£,= J2 b L- (31) 

m— — / 



If the expansion coefficients 6/ m of the temperature fluctuation in the sky are Gaussian, x 2 must obey the x 2 
distribution with 2m+l degrees of freedom. Fig. 9 shows the two non-Gaussian features in the distribution of 6; m 's:a 
slight shift of the peak to the center(zero); slow convergence to zero for large x 2 . As shown in Fig. 10, the distribution 
of x 2 is approximately obtained by assuming that 6; m 's obey G(Z, 1) (actually, the distribution functions of fe; m 's 
are slightly much similar to the Gaussian distributions on large angular scales). The two non-Gaussian features are 
attributed to the nature of the distribution functions of each 6; m which give large values at 6; m ~ and decrease slowly 
at bi m » 1 compared with the Gaussian distributions. 

The slow decrease of the distribution of \ 2 is important in discriminating the non-Gaussian models with the Gaussian 
models. As shown in Fig.ll, observing x 2 ~ 50 are not improbable for the Weeks f^o model (I — 15) whereas it is 
almost unlikely for the Gaussian model. Because the distribution is slowly decreased for large x, the cosmic variances 
(AC;) 2 are expected to be larger than that of the Gaussian models. From Fig. 12, on large angular scales(2 < I < 15), 
one can see that the standard deviations AC; of C; in the two CH models are approximately 1-2 times of that for the 
Gaussian models. 
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FIG. 9. The distributions of x 2 = (2l + l)C'i/Ci for the Weeks model with fi = 0.2, Z = 5(left) and 15(right). The horizontal 
axes represent the values of x 2 ■ The distributions are calculated using 33 eigenmodes (k < 13) based on 200 realizations of the 
initial Gaussian fluctuations <Ev(0), and 200 realizations of the base points. The contribution of modes fc>13 is approximately 
less than 8 percent for ?<15. The solid curves represent the \ 2 distributions with ll(left) and 31(right) degrees of freedom. 




FIG. 10. The distributions of x 2 = (21 + l)Ci/C t in an approximated model in which b; m 's obey G(Z, 1) for 1 = 5 and 1 = 15 
based on 40000 realizations for each 6; m . The horizontal axes represent the values of x 2 - The solid curves correspond to the 
X 2 distributions with ll(left) and 31(right) degrees of freedom. 
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FIG. 11. Plots of 1 — P(Z) (P(Z) is the cumulative distribution function) which gives the probability of observing X > Z. 
The solid curves correspond to 1 — P{x 2 ) f° r the Weeks model Qo = 0.2, I = 5 (left) and I = 15 (right). The dashed curves 
correspond to 1 — P{x 2 ) of the Gaussian model. 
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FIG. 12. Plots of AC; (CH)/ AC; (Gauss) for the two CH models based on 200 realizations of the initial perturbation $„(0) 
and 200 realizations of the base point. ACi denotes the standard deviation (la) of Ci. 



IV. TOPOLOGICAL QUANTITIES 



Topological measures:total area of the excursion regions, total length and the genus of the isotemperature contours 
have been used for testing Gaussianity of the temper ature fl uctuations in the COBE DMR data fl26|,E27|| . Let us first 
summarize the known results for Gaussian fields (see 



29|) 



The genus G of the excursion set for a random temperature field on a connected and simply-connected 2-surface 
can be loosely defined as 



G = number of isolated high-temperature connected regions 
— number of isolated low-temperature connected regions. 



(32) 



For instance, for a certain threshold, a hot spot will contribute +1 and a cold spot will contribute —1 to the genus. 
If a hot spot contains a cold spot, the total contribution to the genus is zero. The genus which is the global property 
of the random field can be related to the integration of the local properties of the field. From the Gauss-Bonnet 
theorem, the genus of a closed curve C being the boundary of a simply-connected region £lc which consists of N arcs 
with exterior angles a\ 1 012, ...oln can be written in terms of the geodesic curvature n s and the Gaussian curvature K 
as 



G = 



2tt 



/ Kgds + V] oti + / KdA 

Jc i=1 Jn c 



(33) 



For a random field on the 2-dimensional Euclidean space E 2 where the N arcs are all geodesic segments (straight line 
segments), K and n g vanish. Therefore, the genus is written as 



1 N 

— Y 



Ge* — 7— > Oti 

i=l 



(34) 



The above formula is applicable to the locally flat spaces such as E 1 x S 1 and T 2 which have E 2 as the universal 
covering space since K and n g also vanish in these spaces. In these multiply-connected spaces, the naive definition 
Eq.(|33"|) is not correct for excursion regions surrounded by a loop which cannot be contracted to a point. 

In order to compute the genus for a random field on a sphere S 2 with radius equal to 1, it is convenient to use a 
map tp:S 2 — {pi} — {p2}^S 1 x (0, n) defined as 

ip : (sin cos 0, sin sin 0, cos#) — > (</>, 6), < 4> < 2vr, < 6 < w, (35) 

where pi and P2 denote the north pole and the south pole, respectively. Because S 1 x (0, it) can be considered as 
locally flat spaces (4>,0) with metric ds 2 — dd 2 + d(f> 2 which have boundaries 8 = 0, ir, the genus for excursion regions 
that do not contain the poles surrounded by straight segments in the locally flat (<j), 6) space is given by Eq. (|34|) . It 
should be noted that the straight segments do not necessarily correspond to the geodesic segments in S 2 . If a pole is 
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inside an excursion region and the pole temperature is above the threshold then the genus is increased by one. If the 
pole temperature is below the threshold, it does not need any correction. Thus the genus for the excursions is 

i 

where 014 is the exterior angles at the intersection of two straight segments in the (</>, 9) space and N p is the number 
of poles above the threshold. 

Now consider an isotropic and homogeneous Gaussian random temperature field on a sphere S 2 with radius 
1. Let (x,y) be the local Cartesian coordinates on S 2 and let the temperature correlation function be C(r) =< 
(AT/T) (AT/T) r >with r = x 2 + y 2 and C = C(0) = a 2 , where a is the standard deviation and C 2 = -(d 2 C/dr 2 ) r=0 . 
Then the expectation value of the genus for a threshold AT/T = va is given as [^9| 

where erfc(:r) is the complementary error function. The first term in Eq.([57|) is equal to the averaged contribution for 
the excursions which do not contain the poles while the second term in Eq.(|37]) is the expectation value of N p . 
The mean contour length per unit area for an isotropic homogeneous Gaussian random field is |2q,[29| 




< s >= - ^ e- v l\ (38) 



and the mean fractional area of excursion regions for the field is the cumulative probability of a threshold level, 

<a>=ier&/-^, (39) 

which gives the second term in Eq. (|37|) . 

As in Sec. Ill, the CMB anisotropy maps for the two CH adiabatic models are produced by using eigenmodes 
k < 13 and angular components 2 < / < 20 for fio = 0.2 and 0.4. The contribution of higher modes are approximately 7 
percent and 10 percent for O = 0.2 and 0.4, respectively. The initial power spectrum is assumed to be the extended 
Harrison-Zel'dovich spectrum. The beam-smoothing effect is not included. For comparison, sky maps for the Einstein- 
de-Sittcr model with the Harrison-Zel'dovich spectrum C; cx 1/(1(1 + 1)) are also simulated. 

In order to compute the genus and the contour length for each model, 10000 CMB sky maps on a 400x200 grid 
in the (</>, 9) space are produced. The contours are approximated by oriented straight segments. The genus comes 
from the sum of the exterior angles at the vertices of the contours and the number of poles at which the temperature 
is above the threshold. The total contour length is approximated by the sum of all the straight segments. Typical 
realizations of the sky map are shown in Fig. 13. 

Fig. 14 and Fig. 15 clearly show that the mean genuses and the mean total contours for the two CH models are well 
approximated by the theoretical values for the Gaussian models. This is a natural result since the distribution of the 
expansion coefficients bi m is very similar to the Gaussian distribution in the modest range. On the other hand, at 
high and low threshold levels, the variances of the total contour lengths and the genuses are larger than that for the 
Gaussian models that can be attributed to the nature of the distribution function of bi m . One can easily notice the 
non-Gaussian signatures from Fig. 16 and Fig. 17. The excess variances for the Weeks model flo = 0.4 compared with the 
Gaussian flat Harrison-Zel'dovich model are observed at the absolute threshold level approximately \v\ > 1.4 for genus 
and \u\ >0.6 for total contour length. If one assumes that the initial fluctuations are given by (<E>„(0)) -2 cx v(v 2 + 1), 
the temperature fluctuations for CH models can be described as Gaussian pseudo-random fields. One can see from 
Fig. 18 that the behavior of the variances of genus and total contour length for the Gaussian CH models is very similar 
to that for the flat Harrison-Zel'dovich model and the variances at high and low threshold levels are considerably 
smaller than that for the non-Gaussian models. 

Because the mean behavior for the two non-Gaussian CH models is well described by the Gaussian models, the 
COBE DMR data which excludes grossly non-Gaussian models cannot constrain the two CH models by the 

topological measurements. However, one should take account of a fact that the signals in the 10° smoothed COBE 
DMR 4-year sky maps are comparable to the noises |3(]] that makes it hard to detect the non-Gaussian signals in 
the background fluctuations. In fact, some recent works using different statistical tools have shown that the COBE 
DMR 4-year sky maps are non-Gaussia n ]3l| - |33| although some authors cast doubts upon the cosmological origin 
of the observed non-Gaussian signals 34 35|. Thus the evidence of Gaussianity in the CMB fluctuations is still not 
conclusive. 
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FIG. 13. Contour maps of the CMB (not smoothed by the DMR beam) for the Thurston model Qo = 0.4 and a flat 
(Einstein-de-Sitter) Harrison-Zel'dovich model Ci oc 1/(1(1 + 1)) in which all multipoles Z > 20 are removed. 
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FIG. 14. The mean genuses averaged over 100 realizations of the initial fluctuations and 100 realizations of the base points 
and ±lcr run-to-run variations at 27 threshold levels for the Weeks and the Thurston models with f2o = 0.2. The dashed curves 
denote the mean values for a Gaussian model where Co and Ci are obtained by assuming that the expansion coefficients of the 
eigenmodes are random Gaussian numbers (the mean is zero and the variance is proportional to v~ 2 \ The solid curves denote 
the mean values for a Gaussian model that are best-fitted to that for CH models. 
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FIG. 15. The mean contour lengths averaged over 100 realizations of the initial fluctuations and 100 realizations of the base 
points and ±lcr run-to-run variations at 27 threshold levels for the Weeks and the Thurston models with f2o = 0.2. The dashed 
curves denote the mean values for a Gaussian model where Co and C2 are obtained by assuming that the expansion coefficients 
of the eigenmodes are random Gaussian numbers (the mean is zero and the variance is proportional to v~ 2 ). The solid curves 
denote the mean values for a Gaussian model that are best-fitted to that for CH models. 
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FIG. 16. The mean genuses and ±la run-to-run variations at 27 threshold levels for a Weeks model with f2o = 0.4 averaged 
over 100 realizations of the initial fluctuations and 100 realizations of the base points and that for a flat Harisson-Zel'dovich 
model averaged over 10000 realizations. The dashed curves denote the mean genuses for the corresponding Gaussian models. 
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FIG. 17. The mean total contour lengths and ±1(7 run-to-run variations at 27 threshold levels for a Weeks model with 
fio = 0.4 averaged over 100 realizations of the initial fluctuations and 100 realizations of the base points and that for a flat 
Harisson-Zel'dovich model averaged over 10000 realizations. The dashed curves denote the mean total contour lengths for the 
corresponding Gaussian models. 
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FIG. 18. The mean total contour lengths and genuses and zbla run-to-run variations at 27 threshold levels for a Weeks model 
with fio = 0.4 averaged over 300 realizations of the base points. Here it is assumed that the initial fluctuations deterministically 
satisfy (<&„(0)) -2 ccv(u 2 + 1) so that the fluctuations are described by the Gaussian statistics. The dashed curves denote the 
mean total contour lengths and the mean genuses for the corresponding Gaussian models. 
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V. CONCLUSION 



In this paper, Gaussianity of the eigenmodes and non-Gaussianity in the CMB temperature fluctuations in two 
smallest CH(Weeks and Thurston) models are investigated. As shown in Sec. II, it is numerically shown that the 
expansion coefficients of the two CH spaces behave as if they are random Gaussian numbers at almost all the places. 
If one recognizes the Laplacian as the Hamiltonian of a free particle, each eigenmode is interpreted as a wavefunction 
in a stationary state. The observed behavior is consistent with a prediction of RMT which has been considered to 
be a good empirical theory that describe the statistical properties of quantum mechanical systems whose classical 
counterparts are strongly chaotic. However, as we have seen, the global symmetries in the system can veil the generic 
properties. For instance, some eigenmodes on the Thurston manifold have a Z2 symmetry at a point where the 
injectivity radius is maximal. For these eigenmodes, the expansion coefficients are strongly correlated;hence they can 
no longer considered to be random Gaussian numbers. 

Because the eigenmodes actually satisfy the periodic boundary conditions, there are points on a sphere S 2 which are 
identified with different points on S 2 . These points form pairs of circles which are identified by the periodic boundary 
conditions |3(J. If one could identify all the circles on a sphere, one would be able to construct the corresponding CH 
space [[37|. Similarly, if one could identify all the fixed points and the corresponding symmetries, one would be able 
to construct a CH space which have these symmetries. The observed "randomness" in the eigenmodes is actually 
determined by these simple structures. 

In order to understand the symmetric structures of the CH spaces, it is useful to choose an observing point (base 
point) at which one enjoys symmetries as many as possible. However, in reality, there is no natural reason to consider 
fluctuations at only these particular points since the CH spaces are globally inhomogeneous. 

Since the CMB fluctuations can be written in terms of a linear combination of eigenmodes, the fluctuations in 
CH models are almost spatially "isotropic" if averaged all over the space except for very limited places at which the 
eigenmodes have certain symmetries provided that the eigenmodes are Gaussian. The spatial "isotropy" implies that 
the contribution of non-diagonal terms in the two-point correlation functions are negligible. Thus the validity of the 
statistical tests using the angular power spectrum Ci [fz| — |l0|| cannot be questioned on the ground that the background 
space is anisotropic at a certain point. 

If one assumes that the initial fluctuations are Gaussian as in the standard inflationary scenarios, the temperature 
fluctuations are described by isotropic non-Gaussian random fields since they are written in terms of a sum of products 
of two independent random Gaussian variables, namely the initial perturbations and the expansion coefficients of the 
eigenmodes. The distribution functions of the expansion coefficients 6/ m for the sky maps at large values are slowly 
converged to zero than the Gaussian distribution with the same variance and the cosmic variances are found to be 
larger than that of the Gaussian models. 

The increase in the variances are much conspicuous for topological quantities at large or small threshold levels. On 
the other hand, the mean behavior is well approximated by the Gaussian predictions. Therefore, the obtained results 
agree with the COBE DMR 4- year maps analyzed in |2(||27j . In real observations one has to tackle with what obscure 
the real signals such as pixel noises, galactic contaminations, beam-smoothing effect and systematic calibration errors 
which have not been considered in this paper. The absence of large deviations from the mean values at large or small 
threshold levels in the current data may be due to these effects, which will be much explored in the future work. 

Although the recent observations seem to prefer the flat FRW models with the cosmological constant, the evidence 
is not perfectly conclusive. If one includes the cosmological constant for a fixed curvature radius, the radius of the last 
scattering surface (horizon) at present in unit of curvature radius becomes large. Therefore the observable imprints of 
the non-trivial topology of the background space become much prominent. For instance, the number Nf of copies of 
the fundamental domains inside the last scattering at the present slice is approximately 27.9 for a Weeks model with 
S1a = 0.6 and Vl m — 0.2 whereas Nf = 4.3 if £!a = and Vt m = 0.& Thus we have still great possibilities in detecting the 
non-trivial topology by the future satellite missions such as MAP and PLANCK which will provide us much better 
information on the statistical properties of the real signals. The large deviations of the topological quantities from 
the mean values would be the good signals that indicate the hyperbolicity (negative curvature) and the finiteness 
(smallness) of the universe in addition to the direct observation of the periodic structures peculiar to each non-trivial 
topology (see |58| for recent developments). 
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